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After an introduction on phenomena due to spin and mass-energy currents on clocks 
and photons, we review the 1995-2001 measurements of gravitomagnetic field of 
Earth and Lense-Thirring effect obtained by analyzing the orbits of the two laser- 
ranged satellites LAGEOS and LAGEOS II; this method has provided a direct mea- 
surement of Earth's gravitomagnetism with accuracy of the order of 20 %. A future 
accurate measurement of the Lense-Thirring effect, at the level of 1% accuracy, may 
include the LARES experiment that will also provide other basic tests of general 
relativity and gravitation. Finally, we report the latest measurement of the Lense- 
Thirring effect, obtained in 2002 with the LAGEOS satellites over nearly 8 years 
of data. This 2002 result fully confirms and improves our previous measurements 
of the Earth frame-dragging: the Lense-Thirring effect exists and its experimental 
value is within ~ 20% of what is predicted by Einstein's theory of general relativity. 
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1 Gravitomagnetic phenomena on test gyroscopes, test particles, clocks 
and photons 

In this paper we review some phenomena arising in the vicinities of a rotating 
body and some proposals and recent measurements of frame-dragging and Lense- 
Thirring effect obtained by laser-ranged satellites. For a general review of tests and 
measurements of general relativistic and gravitational effects we refer to [jl], 0. 

Einstein's general theory of relativity 0] predicts the occurrence of peculiar 
phenomena on test gyroscopes, test particles, clocks and photons in the vicinity of a 
mass-energy current and thus in the vicinity of a spinning body due to its rotation. 

In section (2) we describe how the orbit of a test-particle is influenced by 
the spin of a central body and the direct measurement of the orbital perturbations 
of laser ranged satellites due to the Earth spin, i.e. the Lense-Thirring effect 0. In 
1995-2002 the Lense-Thirring effect was measured with about 20% accuracy using 
the LAGEOS and LAGEOS II satellites § |, |, 0, see section (2). 

Small test-gyroscopes, that determine the axes of a local, freely falling, 
inertial frame, where "locally" the gravitational field is "unobservable" , rotate with 
respect to "distant stars" due to the spin of a body. This effect is described in ^ and 
should be measured with about 1% accuracy by the Gravity Probe-B experiment 
i- 

However, not only test particles and gyroscopes are affected by the spin 
of the central object but also photons and clocks. In this section we review some 
gravitomagnetic phenomena on clocks and photons. A photon co-rotating around a 
spinning body takes less time to return to a "fixed point" (with respect to distant 
stars) than a photon rotating in the opposite direction |1|, 0. Since light rays are used 
to synchronize clocks, the different travel-time of co-rotating and counter-rotating 
photons implies the impossibility of synchronization of clocks all around a closed 
path around a spinning body; the behavior of light rays and the behavior of clocks 
around a spinning body are intimately connected. In several papers the "frame- 
dragging clock effect" around a spinning body has been estimated and some space 
experiments have been proposed to test it jT^, |rT|, ^. Thus, when a clock, co-rotating 
very slowly (using rockets) around a spinning body and at a constant distance from 
it, returns to its starting point, it finds itself advanced relative to a clock kept there 
at "rest" (with respect to "distant stars", see above). Similarly a clock, counter- 
rotating arbitrarily slowly and at a constant distance around the spinning body, finds 
itself retarded relative to the clock at rest at its starting point ^ |[. For example, 
when a clock that co-rotates very slowly around the spinning Earth, at r ~ 6000 
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km altitude, returns to its starting point, it finds itself advanced relative to a clock 
kept there at "rest" (with respect to "distant stars") by At ~ ^^^^ ~ 5 x 10^"'^'^ 
sec, where = 145 cm^ is the Earth angular momentum. Similarly, a clock, that 
counter-rotates very slowly around the spinning Earth, finds itself retarded relative 
to a clock kept there at "rest" by the same amount. Then, the difference between 
the time read by the two clocks when they meet again after a whole revolution is 
about ~ 10-16 [g g. 

However, Einstein's gravitational theory predicts peculiar phenomena also 
inside a rotating shell |]l|. In reference ^ we derive the time-delay in travel-time 
of photons due to the spin of a body both outside a rotating body and inside a 
rotating shell. We then show that this time-delay by the spin of an astrophysical 
object might be detected in different images of the same source by gravitational 
lensing. 

Since here we are only interested to analyze the time delay due to spin, we 
chose a simple configuration where source, lens and observer are aligned and we use 
quasi-Cartesian coordinates 0. We then get, for a photon with impact parameter 
b traveling on the equatorial plane of the source: AT = 4Mln (^2|j -|- In this 
expression ^ ^ 6 is the distance of source and observer from the lens, the first 
term is the standard Shapiro time delay and the second term is the gravitomagnetic 
time-delay due to the spin of the deflecting body. 

Let us now give the order of magnitude of the time delay due to the spin 
of some astrophysical sources. For the sun, by considering two light rays on the 
equatorial plane of the Sun, grazing the Sun on opposite sides, the relative grav- 
itomagnetic time delay is AT/^i = 3.35 ■ 10"^^ sec. For the lensing galaxy of the 
Einstein cross, by assuming a simple model for rotation and shape of the central 
object (see and references therein), we then get: AT{^ ~ 8 hr. As a third example 
we consider the relative time delay of photons due to the spin of a typical cluster 
of galaxies. Depending on the geometry of the system and on the path followed 
by the photons, we then find relative time delays ranging from a few minutes to 
several days. Then, at least in principle, one could detect the time delay due to the 
spin of a lensing galaxy by removing the larger quadrupole-moment time delay by 
a method described in of course, as in the case of the Sun, one should be able 
to accurately enough model and remove all the other delays, due to other physical 
effects, from the observed time delays between the images. 

Let us now analyze the time delay in the travel time of photons propagating 
inside a shell of mass M rotating with angular velocity u. Inside a spinning shell it 
is in general not possible to synchronize clocks all around a closed path. Indeed, if 
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we consider a clock co-rotating very slowly on the equatorial plane along a circular 
path with radius r, when back to its starting point it is advanced with respect to 
a clock kept there at rest (in respect to distant stars). The difference between the 
time read by the co-rotating clock and the clock at rest is equal to: 5T = j^nur"^. 
For a shell with finite thickness we just integrate this expression from the smaller 
radius to the larger one. 

If we now consider a photon traveling with an impact parameter r on the 
equatorial plane of a galaxy; the time delay due to the rotation of the external mass 
for every infinitesimal shell with mass dm = ATipR'^dR' and radius R' > \r\, is 0: 
Atdm = ^^ujr ^'^~^^ . This is the time delay due to the spin of the external thin 
shell. By integrating this expression from |r| to the external shell radius R, we have: 

AT = —cur r pR'VR'^ -r^dR' (1) 

3 J\r\ 

This is the time delay due to the spin of the whole rotating mass of the 
external shell. From this formula we can easily calculate the relative time delay 
between two photons traveling on the equatorial plane of a rotating shell, with 
impact parameters ri and r2. 

Let us give the time delay corresponding to some astrophysical configu- 
rations. In the case of the "Einstein cross", after some calculations, based on a 
standard model for the lensing galaxy (see ref. and references therein), the order 
of magnitude of the relative time delay of two photons traveling at a distance of 
ri ~ 650 pc and r2 ^ — 650 pc from the center, using (|l]) in the case ri ~ — r2, is: 
AT ~ 20 min. If the lensing galaxy is inside a rotating cluster, or super-cluster, to 
get an order of magnitude of the time delay due to the spin of the mass rotating 
around the deflecting galaxy, we use typical super-cluster parameters (see and 
references therein). If the galaxy is in the center of the cluster and light rays have 
impact parameters ri ~ 15 kpc and r2 — — 15 kpc (of the order of the Milky Way 
radius), the time delay, applying formula (|1|) in the case ri ~ — r2 and p = constant, 
is: At ~ 1 day. 

Finally, if the lensing galaxy is not in the center of the cluster but at a 
distance r = aR from the center, with < a < 1 and R radius of the cluster, by 
integrating (|l]) between r = aR and R, when ri ~ r2 we have: 

AT = ^coir^ - r2)p(l - a^Y^^l - Aa')R' (2) 
9 

Thus, if the lensing galaxy is at a distance of 10 Mpc from the center of 
the cluster, the relative time delay due to the spin of the external rotating mass 
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between two photons with (ri — — 30 kpc, is: AT ~ 0.9 day. Since the present 



measurement uncertainty in the lensing time delay is of the order of 0.5 day |]T2 
the spin time delay might already be observable. 

In conclusion, in ref. we have derived and studied the " spin-time- delay" 
in the travel time of photons propagating near a rotating body, or inside a rotating 
shell due to the angular momentum. We found that there may be an appreciable 
time delay due to the spin of the body, or shell, thus spin-time- delay must be taken 
into account in the modeling of relative time delays of the images of a source ob- 
served at a far point by gravitational lensing. This effect is due to the propagation 
of the photons in opposite directions with respect to the direction of the spin of the 
body, or shell. If other time-delays can be accurately enough modeled and removed 
from the observations 0, one could directly measure the spin-time-delay due to the 
gravitomagnetic field of the lensing body. We have analyzed the relative time-delay 
in the gravitational lensing images caused by a typical rotating galaxy, or cluster 
of galaxies. We have then analyzed the relative spin-time-delay when the path of 
photons is inside a galaxy, a cluster, or super-cluster of galaxies rotating around the 
deflecting body; this effect should be large enough to be detected at Earth. The 
measurement of the spin-time-delay, due to the angular momentum of the external 
massive rotating shell, might be a further observable for the determination of the 
total mass-energy of the external body, i.e. of the dark matter of galaxies, clusters 
and super-clusters of galaxies. Indeed, by measuring the spin-time-delay one can 
determine the total angular momentum of the rotating body and thus, by estimating 
the contribution of the visible part, one can determine its dark- matter contribution. 
The estimates presented in this paper are preliminary because we need to apply the 
spin-time-delay to some particular, known, gravitational-lensing images. Further- 
more, we need to estimate the size and the possibility of modeling other sources 
of time-delay. Nevertheless, we conclude that, depending on the geometry of the 
astrophysical system considered, the relative spin-time-delay can be a quite large 
effect. 



2 Measurement of gravitomagnetism with laser ranged satellites 

In ref. (see also ref. |T^) we describe the LARES experiment to measure 

the Lense-Thirring effect with relative accuracy of about 1%. This laser ranged 
satellite, by detecting its perigee rate, would also test the foundations and other basic 
phenomena of general relativity and gravitational interaction. Indeed, LARES would 
improve the bounds on hypothetical long-range gravitational forces and the bounds 
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on deviations from the inverse square law for very weak-field gravity; LARES would 
improve, by about two orders of magnitude, the accuracy in testing the equivalence 
principle and would provide an improved measurement in the field of Earth of the 
PPN (Parametrized- Post-Newtonian) parameters ai, (3 and 7 |]13|. 

In this section we describe the 1995-2002 measurements of the Lense- 
Thirring effect obtained using LAGEOS and LAGEOS II. 

The measurement of distances has always been a fundamental issue in as- 
tronomy, engineering, and science in general. So far, laser ranging has been the most 
accurate technique for measuring the distances to the moon and artificial satellites. 
Short-duration laser pulses are emitted from lasers on Earth, aimed at the target 
through a telescope, and then reflected back by optical cube-corner retrorefiectors 
on the moon or an artificial satellite |15|, such as LAGEOS. By measuring the total 
round-trip travel time, one can determine the distance to a retroreflector on the 
moon with an accuracy of about 2 cm and to the LAGEOS satellites with a few 
millimeters accuracy. 

Our detection and measurement of the Lense-Thirring effect was obtained 
by using the satellite laser ranging data of the satellites LAGEOS (LAser GEO- 
dynamics Satellite, NASA) and LAGEOS II (NASA and ASI, the Italian Space 
Agency) and the Earth gravitational field models, JGM-3 and EGM-96. The LA- 
GEOS satellites are heavy brass and aluminum satellites, of about 406 kg weight, 
completely passive and covered with retrorefiectors, orbiting at an altitude of about 
6,000 km above the surface of Earth. LAGEOS, launched in 1976 by NASA, and 
LAGEOS II, launched by NASA and ASI in 1992, have an essentially identical struc- 
ture but they have different orbits. The semimajor axis of LAGEOS is a = 12, 270 
km, the period P = 3.758 hr, the eccentricity e = 0.004, and the inclination 
/ = 109.9°. The semimajor axis of LAGEOS II is ajj = 12, 163 km, the eccentric- 
ity eii = 0.014, and the inclination /// = 52.65°. 

We analyzed the laser-ranging data adopting the lERS (International Earth 
Rotation Service) conventions [16| in our modeling, however, in the 1998 analysis, we 
used the static and tidal EGM-96 model [0. Error analysis of the LAGEOS orbits 
indicated that the EGM-96 errors can only contribute periodic root-sum-square er- 
rors of 2 to 4 mm radially, and in all three directions they do not exceed 10 to 17 mm. 
The initial positions and velocities of the LAGEOS satellites were adjusted for each 
15-day batch of data, along with variations in their reflectivities. Solar radiation 
pressure. Earth albedo, and anisotropic thermal effects were also modeled [|T^]. In 
modeling the thermal effects, the orientation of the satellite spin axis was obtained 
from ref. M]. Lunar, solar, and planetary perturbations were also included in the 
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equations of motion, formulated according to Einstein's general theory of relativity 
with the exception of the Lense-Thirring effect, which was purposely set to zero. 
All of the tracking station coordinates were adjusted (accounting for tectonic mo- 
tions) except for those defining the TRF terrestrial reference frame. Polar motion 
was also adjusted, and Earth's rotation was modeled from the very long baseline 
interferometry-based series SPACE96 PD[. We analyzed the orbits of the LAGEOS 
satellites using the orbital analysis and data reduction software GEODYN II |ETI . 



The node and perigee of LAGEOS and LAGEOS II are dragged by the 
Earth's angular momentum. From the Lense-Thirring formula 0, ^, P], we get 

r\Lense—Thirrinq rvj 01 / J ^, Lense—Thirring ^ q-, r 1 . .11. 1 

ill = ol mas/yr and iln = ol.o, where mas is a miUisecond 

of arc. The argument of pericenter (perigee in our analysis), uj, also has a Lense- 
Thirring drag thus, we get for LAGEOS: ^^^"*^-™*'^''*"'f ^ 32 mas/yr, and for 
LAGEOS II: (^^^'^■'^-™*''™9 ^ — 57 mas/yr 0, |^, |^. The nodal precessions of 
LAGEOS and LAGEOS II can be determined with an accuracy of the order of 1 
mas/year. Over our total observational period of about 4 years, we obtained a root 
mean square (RMS) of the node residuals of about 4 mas for LAGEOS and about 
7 mas for LAGEOS II 0. For the perigee, the observable quantity is the product 
e-a-cj, where e is the orbital eccentricity of the satellite. Thus, the perigee precession 
Lo for LAGEOS is difficult to measure because its orbital eccentricity e is ~ 4 x 10~^. 
The orbit of LAGEOS II is more eccentric, with e ~ 0.014, and the Lense-Thirring 
drag of the perigee of LAGEOS II is almost twice as large in magnitude as that of 
LAGEOS. Over about 4 years, we obtained a root mean square of the residuals of 
the LAGEOS II perigee of about 25 mas [Q, whereas the total Lense-Thirring effect 
on the perigee is, over 4 years, = —228 mas. 

To precisely quantify and measure the gravitomagnetic effects we have 
introduced the parameter /i that is by definition 1 in general relativity and zero 
in Newtonian theory. 

The main error in this measurement is due to the uncertainties in the 
Earth's even zonal harmonics and their time variations. The unmodeled orbital 
effects due to the harmonics of lower order are comparable to, or larger than, the 
Lense-Thirring effect. However, by analyzing both the JGM-3 and the EGM-96 
models with their uncertainties in the even zonal harmonic coefficients and by cal- 
culating the secular effects of these uncertainties on the orbital elements of LAGEOS 
and LAGEOS II, we find 0] that the main sources of error in the determination of 
the Lense-Thirring effect are concentrated in the first two even zonal harmonics, 
J2 and J4. We can, however, use the three observable quantities fi/, Vtn and ton 
to determine jU [Q, thereby avoiding the two largest sources of error, those arising 
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from the uncertainties in J2 and J4. We do this by solving the system of the three 
equations for the three residuals 517/, dVLn and duu in the three unknowns /i, J2 
and J4, obtaining: 

^ ^LAGEOSI + Cl 5 ^LAGEOSII + ^2 ^ ^LAGEOSII — 



= /i(31 + 31.5 Ci — 57 C2)mas/yr + other errors = fi {60.2 mas/ yr), (3) 

where Ci = 0.295 and C2 = —0.35. Equation (^) for /i does not depend 
on J2 and J4 nor on their uncertainties; thus, the value of fi that we obtain is 
unaffected by the largest errors, which are due to 6J2 and 6J4, and is sensitive only 
to the smaller errors due to SJ2n, with 2n > 6. 

Similarly, regarding tidal, secular, and seasonal changes in the geopotential 
coefficients, the main effects on the nodes and perigee of LAGEOS and LAGEOS II, 
caused by tidal and other time variations in Earth's gravitational field [H, |^, |14|, are 



due to changes in J2 and J4. However, the tidal errors in J2 and J4 and the errors 
resulting from other unmodeled, medium and long period, time variations in J2 and 
J4, including their secular and seasonal variations, are eliminated by our combination 
of residuals of nodes and perigee. In particular, most of the errors resulting from 
the 18.6- and 9.3-year tides, associated with the lunar node, are eliminated in our 
measurement. An extensive discussion of the various error sources that can affect 
our result is given in 0, ^ |l^, only a brief discussion of the error sources is given 
below. 

Let us now report the main results of our measurements. In Fig. 1, we 
display the improved analysis (obtained with the linear combination of the resid- 
uals of the nodes of LAGEOS and LAGEOS II and perigee of LAGEOS II according 
to Eq. (0)) using the static and tidal Earth gravitational model EGM-96, we also 
refined the non-gravitational perturbations model, the total period of observations 
was of 4 years, longer of about 1 year with respect to the previous observational pe- 
riod corresponding to ref. We only removed three small periodic residual signals 
and the small observed inclination residuals. The removal of the periodic terms was 
achieved by a least squares fit of the residuals using a secular trend and three pe- 
riodic signals with 1044-, 820-, and 569-day periods, corresponding to, respectively, 
the nodal period of LAGEOS, and the perigee and nodal periods of LAGEOS II. The 
820-day period is the period of the main odd zonal harmonics perturbations of the 
LAGEOS II perigee; the 1044- and 569-day periods are the periods of the main tidal 
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orbital perturbations, with 1=2 and m=l, which were not ehminated using Eq. (^). 
Some combinations of these frequencies correspond to the main non-gravitational 
perturbations of the LAGEOS II perigee. We notice that this analysis, using EGM- 
96 and its accurate tidal model, is substantially independent of the removed signals, 
whereas the previous analysis 0, was in part sensitive on the periodic terms in- 
cluded in the fit. In other words, our value (Fig. 1) of the secular trend does not 
significantly change by fitting additional periodic perturbations, and indeed, even 
the fit of the residuals with a secular trend only, with no periodic terms, just changes 
the slope by about 10 %. Nevertheless, in this case, the root mean square of the 
post-fit residuals increases by about four times with respect to Fig. 1. 

milliarcsec 

300 




seconds 



Figure 1: 1997-measurement of the Lense-Thirring effect. Combination of the resid- 
uals of the nodes of LAGEOS and LAGEOS II and perigee of LAGEOS II according 
to Eq. 1^, using the Earth gravitational model EGM-96, over a 4-yGar period j^. 

Our best-fit straight lines of Fig. 1, through the combined residuals of 
nodes and perigee, has a slope ^^easwred ^ -^ -^ _[_ q_q3^ where 0.03 is the standard 
deviation of the fit. This combined, measured, gravitomagnetic perturbation of the 
satellites' orbits corresponds, in a 4-year period, to about 16 m at the LAGEOS 
altitude, that is, about 265 mas. 

The root mean square of the post-fit combined residuals corresponding to 
Fig. 1 is about 9 mas. Our total systematic error is estimated to be of the order of 
30%-40% of ^cB. corresponding to the previous analyses of ref. 0, and of the order 
of 20%-25% of fiGR corresponding to Figure 1 0. 

Using the JGM-3 covariance matrix, we found the errors due to the un- 
certainties in the even zonal harmonics J2n, with 2n > 6, to be: ^^e''en2onaZs:J2„>J6 
~17% of fiGR, and using the EGM-96 covariance matrix: 5^e''en2onais:j2„>j6 <i3% 
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of jJ,GR- The errors in the modehng of the perigee rate of LAGEOS II due to the un- 
certainties in the odd zonal harmonics J2n+i are, with EGM-96: 6 ^°'^^°'>^'^'-^ ^2 % of 
fJ^GR- Using the EGM-96 tidal model, we estimated the effect of tidal perturbations 
and other variations of Earth gravitational field to be S ^t'^des + othervanations 
^GR- On the basis of analyses |22| of the non-gravitational perturbations — in 
particular, those on the perigee of LAGEOS II — we found 5 ^rion- gravitational < 13%. 
20% of /^Gi?! including uncertainties in the modeling of the satellites' reflectivities, 
and the error due to uncertainties in the orbital inclinations of LAGEOS and LA- 
GEOS II was estimated to be 5 ^iricUnation <^ % 

Taking into account all of these error sources, we arrived at a total root- 



sum-square error ~2U7o-25% of Hgr- Therefore, over an observational period of 4 
years and using EGM-96, we determined ^A^e"*"''^'^ = x_x -|- o.25 0, where 0.25 is 
the estimated total uncertainty due to all the error sources — see Fig. 1. 
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2002-Measijrement of the Lense-Thirring Effect 



Figure 2: Latest, 2002, measurement of the Lense-Thirring effect using LAGEOS 
and LAGEOS II, obtained by only modeling the radiation pressure coefficient of 
LAGEOS II, over nearly 8 years of data. The best-fit line shown through these 
combined residuals has a slope /i ~ 1 ± 0.02. The total estimated systematic error 
is about ± 0.2. The total measured signal is nearly 440 milliarcsec and the RMS of 
the post-fit residuals is about 10 milliarcsec 
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We finally briefly report on our latest, 2002, measurement of the Lense- 
Thirring effect over 7.3 years of data of LAGEOS and LAGEOS 11, i.e. over an 
observational time nearly double than the longer period of our previous analyses, 
obtained by only modeling the radiation pressure coefficient of LAGEOS 11 (see Fig. 
2) i- 

This recent measurement fully confirms and improves our previous results: 
the Lense-Thirring effect exists and its experimental value, /i = (1±0.02)±0.2 (where 
±0.02 is the standard deviation of the fit and ±0.2 is the estimated total systematic 
error), fully agrees with the general relativity prediction of frame-dragging. It is 
important to notice that: (1) in the analysis corresponding to Fig. 2 we only modeled 
on LAGEOS 11 the radiation pressure coefficient of the satellite, i.e. the refiectivity 
coefficient, Cr, and no other parameters such as the along-track accelerations as in 
our previous analyses corresponding to Fig. 1 and ref. [§, P); (2) the RMS of the 
residuals corresponding to Fig. 2 is about 10 milliarcsec whereas the total measured 
signal is about 440 milliarcsec, and finally (3) the quality of the fit and of the 
corresponding measurement can be further improved by further reducing the RMS 
of the 15-day fits (corresponding to each point of Fig. 2) with further processing of 
the data using GEODYN/SOLVE, thus further reducing the RMS of the final fit of 
Fig. 2. 

In conclusion, the Lense-Thirring effect exists and its experimental value, 
/i = 1 ± 0.2, fully agrees with the prediction of general relativity 
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